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WEBER’S THEORY OF THE KERNLEITER 


ALVIN M. WEINBERG 
THE UNIVERSITY OF CHICAGO 


The potential distribution about a kernleiter is determined ac- 
cording to Weber’s method. It is shown that the distribution re- 
duces to the solution of a telegrapher’s equation when the volume of 
the external medium is small. The velocity of propagation as a func- 
tion of the external volume is determined approximately. This in- 
volves the solution of the equation 


Lyhh Spee lees (2)1. 
[J,(ké)]’  [&-4 J, (€)]’ 


where Y, and J, are Bessel functions, and roots of this equation are 


tabulated. The velocities thus found reduce to Lillie’s values deter- 
mined experimentally on the iron wire when the conducting medium 
liso Deviations from these values are predicted for larger 
volumes. 


Most of the mathematical discussions of the kernleiter* have be- 
gun with the assumption that the kernleiter could be replaced by an 
electrical network in which the transverse and longitudinal imped- 
ances had certain known values per centimeter. (See, however, Crem- 
er, 1899; Bishop, 1937; Wilson, Macleod, and Barker, 1933). Such a 
network could be described by a telegrapher’s equation, the solution 
of which gives the current and voltage distribution along the line. 

This classical solution of the cable problem is one-dimensional, 
since it tells nothing about the radial distribution of currents at dis- 
tances out from the kernleiter. Such information becomes important, 
however, in discussing the theory of one of R. S. Lillie’s passive iron 
wire experiments (1925, 1936). Lillie enclosed a passive iron wire 
in a glass cylinder of radius a; he then determined the velocity of pro- 
pagation as a function of a. Since the action currents spread radially 
throughout the space between the wire and its enclosing jacket, one 
would suppose that the velocity of propagation should depend on the 
nature of this radial distribution. Lillie, in analyzing his experimental 
results, assumed that the current density was uniform throughout the 
space between the wire and the jacket. This is correct if, as in his 

* By a “kernleiter” (core-conductor) we mean a cylindrical conductor cov- 


ered by a sheath or membrane and immersed in an electrolytic medium. Nerve 
axon, the passive iron wire, and submarine cable fall into this category. 
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experiments, the jacket is not much larger than the wire; but if the 
radius a is much greater than the wire radius, the current density will 
not be uniform, but will fall off radially from the wire. 

In this paper we shall determine the radial distribution of the 
currents, and we shall then calculate the velocity of an impulse inso- 
far as it depends on this distribution. In our discussion, we shall fol- 
low the treatment of H. Weber (1873a) rather closely. Weber’s work 
has long been recognized as a standard contribution by mathemati- 
cians and physicists (Watson, p. 787), but it seems to have been some- 
what neglected by physiologists, perhaps because, at the time it ap- 
peared, the refined experimental techniques required to check his theo- 
retical predictions were not available. 


Stationary Distribution: Finite External Medium 


We shall first find the stationary current and voltage distribution 
about a semi-infinite kernleiter of radius b covered with a sheath or 
membrane of effective conductivity 1/7 in ohm cm, and surrounded 
by a glass cylinder of radius a, the intervening medium being filled 
with an electrolyte of specific resistance x (Figure 1). For the passive 


ELECTROLYTE 
RESISTIVITY, 2 / 


\ 
SHEATH RESISTIVITY Z,z’ 
. 


FIGURE 1 


iron wire, the electrolyte, 70% HNO;, has a conductivity of about 1 
ohm- cm-; this is only one millionth the conductivity of iron, and so 
the core of the wire will be assumed to be a perfect conductor at zero 
potential. Also, we shall assume that the current distribution is ra- 
dially symmetric; the more general case in which the distribution de- 
pends on the azimuthal angle is treated in detail by Weber. 

If uw is the potential at some point in the medium, then u must 
satisfy Laplace’s equation in cylindrical coordinates: 


O(r as 
or Cn) a 0 1 
LO a2 (1) 
with the following boundary conditions: 
At r =a, no radial current flow: 
OU 


=") (2) 
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At z= o, potential is zero: 


u=O0. (3) 


At r= b, current through the membrane sheath is equal to the 
current reaching the membrane from the electrolyte: 


where h = 7/x. This boundary condition corresponds to Kirchhoff’s 
current law as used in the derivation of the telegrapher’s equation. 

To complete the specification of the problem, the initial distribu- 
tion of the current or voltage must be given. As an approximation to 
the physical situation that obtains when an annular element on the 
kernleiter surface is excited, we may suppose that at z = 0, the longi- 
tudinal current through the sheath is just the action current, while 
all other currents in the system are zero: i.e., atz=0, 


ce ere ALE, (5) 


where 
®(¢) =lim (7J/22bAb) for b<r<b+Ab; 
Ab-0 


Dry 09 Mort eo Ab: 


I being the action current, A b the thickness of the sheath through 
which the current flows, and 7’ the specific resistance of the sheath in 
ohm cm. This, of course, is only an approximation to the actual initial 
condition; it will be seen presently that the exact form of this assump- 
tion is not very important for certain of our purposes. (The function 
u defined by (1)-(5) is essentially the Green’s function (Smythe, p. 
56) for the enclosed kernleiter.) 
The solution of this problem is, according to Weber, 


u(r, Zz) = A; f(x; r) Cae (6) 


where ( after correcting Weber’s normalizing factor) 
f (25,7) =Jo(%i 7) Y1(x1a) — Ji (4, 2) Y.(%: 7), 


Ap i P(r) f(ui,7r) rdr / Xi J Pai, 7ryrdr (7) 


=e ite Li 7 LiAes, b) /[4h? Se b2(1 + h? Bei (es ; b)] ’ 
and the x; are roots of the characteristic equation 
Y,(ax) hx Y,(bx) + Yo(bx) 


=> (8) 
J,(az) haJd,(bx) + Jo(bx) 
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the functions Ji(bx) and Y;(ba) being Bessel functions of the first 
and second kinds, respectively (Watson, pp. 38, 63). The series (6) 
diverges at z = 0, but converges for z > 0; such behavior is charac- 
teristic of Green’s functions. 

The potential at 7 = b , which is just the voltage across the mem- 
brane, is found from (6) and (7): 


u(b,2) =ah? sr TD a; f2(a,b) eo /4h? 


9 
= Nae b2(1 + h? 4,7) f? (xi, 0). ( ) 
If the x; increase sufficiently rapidly with 7, or if z is large, only the 
first term of the expansion (6) need be kept: 
a h? rel £4 f (41,6) 
y See 2 2 ee ae eee L) E-*% 10 
Us (7 me) Ah? — x2 b2(1 = h2 2) Pdi, b) ) ( ) 
or 
1 (bo 2) — tae (11) 


where physically u, must be approximately the potential on the sur- 
face of the kernleiter at z = 0, the approximation improving as the 
ratio x,/x;, 7 > 1, decreases. Since wu has physical significance inde- 
pendent of any artificial assumptions about ®(7), it is clear that to 
this degree of approximation the original initial condition is imma- 
terial. It is seen from (10) or (11) that the potential falls off ex- 
ponentially with length constant or characteristic length 1/2; this 
result corresponds to the exponential distribution found along a cable 
characterized by a telegrapher’s equation. 


The Characteristic Equation 


The characteristic equation (8) was not discussed at any length 
by Weber, and so we shall indicate in some detail how its roots may 
be approximated. If h is very small, (8) becomes 


Yi(ax)  Y.(bx) 


Ji(ax)  Jo(bx) ’ pe 
while if h is very large, the equation reduces to 
Yi(ax) Y.i(bzx) 
1 _ 1 (13) 


Ji(ax) — F,(bx) ° 


Series for the roots of both these equations have been given in 
the literature (Kalaihne, 1907; Sasaki, 1914), and the first few roots 
of (13) have been tabulated by Kalahne for various values of k(=a/b). 
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A similar series for roots of (8) can be obtained by writing (8) as 


[Y.(ax)]’ [x Y,(bx)]' 
[Jo(ax)]’ [x J, (bx) ]' 


(where primes denote derivatives with respect to x), and applying 
Stokes’ method (McMahon, 1894) ; but the resulting series converges 
too slowly and is not convenient for calculation when k is large. 

A rather neat method used by Kalahne to obtain roots of (13) 
can be modified to find as many roots of (8) as are needed for the 
expansion (9). Upon putting 7 = ax, = bx, we can write (8) 


Yi(y) — R(E, b/h) J, (m) =0 (15) 


(14) 


where 
hE Y,(€) +0 Y,(é) 
WEIN EY OIu(E) 


For any assigned value of €, R(é , b/h) is a function of b/h 
alone, Rz(b/h). Then for every value of b/h, Re(b/h) has some con- 
stant value, Res, , so that (15) becomes 


Y1(y) — Reon Ji(y) =0. (17) 


Equation (17), for a preassigned value of &; , is a fairly simple trans- 
cendental equation whose root of the j-th order, 7; , can be found with 
considerable accuracy from Watson’s tables (where Y, and J, are 
listed in adjoining columns), and with greater accuracy from the 
British Association tables (1937). In this way a root of (8) 


R(é,b/h) = (16) 


x, = §:/b (18) 
is found for a value of k given by 
k=n;/&:. (19) 


Roots for intermediate values of k can be found by interpolation. 
We shall now prove that for a fixed k and b/h, the order of “7; as 
a root of (8) is the same as the order of &; as a root of 


Ron (éi) — Reon =O. (20) 


To do this, we must first show that the values of &; which satisfy (20) 
form a monotonic sequence in which the subscript 7 denotes the period 
(defined as a connected range of & over which Ry,(&) takes every ar- 
bitrary real value once and only once) in which &; is found. These 
periods consecutively cover the entire positive real é-axis, since we can 
show that the poles and zeros of R,,(&) are interlaced (Watson, p. 
479) ; this is true since (cf. equations (24), below) 
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dR 2 fle hé 
Soe eee i a OD PACS 
easel pitt] / [Faw +4) 
being always positive for € > 0, implies that Ry,,(&) is monotonically 
increasing. Therefore, the zeros and poles of R are interlaced, and 
the entire positive — — axis is covered by periods in each of which a 
solution of (20) occurs once and only once (Figure 2). The justifica- 


2 
? 


Rul) 


FIGURE 2 
tion for referring to i as the “order” of &; as a root of (20) lies in this 
result. 
We can now show that the order of x; as a root of (8) is the same 


as the order of &; as a root of (20) if we can prove that for a fixed k, 
the function 


¥(€) = Yi(kK&) — Ron (&) Ja (k &) 


is monotonic in (= #/b). For if Y is monotonic, its poles will again 
divide the € — axis into periods over each of which ¥ monotonically 
takes on all real values. But the poles of Y are identical with the poles 
of R(&), and so the i-th zero of ¥ and the i-th zero of Ron(&) — Ream 
must occur in the same period. 

To prove that ¥(&) is monotonic, we first suppose that ¥(é) has 
a horizontal tangent and then show that this leads to a contradiction. 
If ¥ has a horizontal tangent, then for at least one value of b/h, 
must have a multiple root. But since the functions f(#:,7) and 
f(a;,7) are orthogonal (Weber, 1873a) whenever x; # x; but are 
identical if x; = x; , and since this is true regardless of the value of 
b/h, it follows that ¥ cannot have a multiple root for any value of 
b/h, and so ¥ must be monotonic. 

In actually calculating the roots, we start with an arbitrary small 
value of € and ascertain its order as a root of (20) by inspecting Wat- 
son’s tables. Corresponding values of 7 are found by trial, and from 
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these, values of k can be computed [equation (19) ]. First order roots 
computed in this way for several values of k and b = .1588 cm. (1/16 
inch), h=1cm., are shown in Table I. Higher order roots for k = 20 
and the same values of b and h are given in Table II, and x; versus k 
eee which were used in constructing Table II are given in Figure 


22 


° 2 4 6 8 10 12 14 16 6 20' 22 24 26 28 30 
K=o/p—> 


FIGURE 3 
Characteristic values, ~; versus k. Note that x; — ma, increases as k decreases. 


TABLE I: First order roots of (8) compared with x, 
fe — Hem. 6 —.1585: em: 


k—a/b L %, — [2b/h (a2 — b2)]# 
1.020 17.63 17.59 
1.039 12.59 12.58 
1.147 6.30 6.32 
1.50 3.15 3.18 
1.705 2.52 2.57 
1.91 2.14 2.19 
2.10 1.89 1.94 
2.38 1.64 1.69 
2.88 1.26 1.32 

, 3.93 0.882 .936 
5.30 0.630 681 

12.45 0.252 .286 

23.90 0.126 149 


TABLE II: Higher order roots of (8) for k = 20, 
h=1cm., b = .1588 cm. 


k G4 Lo Le X, Xs Xe Xn tg 
Oe ge ie oo PROP AS bee Gy. 78, 
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If the difference a — b is small (i.e. & ~ 1) the current density 
should be practically uniform throughout the external medium. The 
longitudinal resistance per centimeter in this case is x/a(a? — 0b?) 
and so xa, the reciprocal of the characteristic length as calculated 
from the usual telegrapher’s equation (Weinberg, 1940) is 


Te ene 
ean ii a ct oe a 


where M is the longitudinal resistance per centimeter, and 1/s is the 
transverse conductance per centimeter. The first order root of (8) 
must reduce to xa in this case; we can prove this as follows: 


Ifa=b+4, (& << 1), then by Taylor’s theorem 
Y, (ax) =Y,(bu + Ax) © Y,(b2) + Ax Yy' (ba) ee 1" (bx) 
oe (22) 
J, (aa) =I,(bu + Aa) ©I(bx) + Ax Is (bu) + ane (ba) . 
If these expressions are substituted into (8) we find, after clearing 


fractions and using certain identities between Bessel functions of the 
first and second kinds, 


h A? 1 2A ey TA 
ae A — Zhi a = 
o( 7 L )+ 5 (2 7 ois Bs 0. (23) 


The identities used in obtaining (23) are (Watson, p. 76) 
In (E) V'n(E) — Yn(E) In (E) = 2/n& 
In(€) V'n(S) — Yn(E) In (€) = — 2/n (24) 


TnlEV" ACE) — YaCSI"9() => (“PF 1). 


Equation (23) is a quadratic in « which has the positive root 
eee 2? — b A Hz) 
Ab(2hb—hA+ BA) 


this being a second approximation to the first root of (8). If (A /b)? 
be neglected in comparison with 4/b , the root becomes 


= 20 — @ 
\ h(a — b) (26) 


(25) 
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‘ . atdb 
and since lim ——_ — b , 
A>0 


lim 2 = %a. 

A>0 
In the limit, the length constant in the exact theory is identical with 
that of the simple telegrapher’s theory [cf. Carslaw (1906), p. 322, 
and Table II above]. Moreover( Figure 3), the ratio x,/x;(i > 1) de- 
creases with k so that whenk ~ 1, only the first term of (9) remains, 
and the solution of the present partial differential equation reduces to 
the solution of the simple telegrapher’s equation. 

We should expect the fall of potential along the line to be more 
gradual for a uniform radial current distribution than for the actual 
current distribution, since the rate of fall of potential is measured by 
the effective longitudinal conductance [equation (21)], and the effec- 
tive conductance is overestimated by assuming a uniform distribu- 
tion. Table I seems to contradict this result, for x. [which by equa- 
tion (21) is inversely proportional to the longitudinal conductance] 
is larger than the corresponding x,. The discrepancy becomes worse 
as k increases; that is, as the non-uniformity of the radial distribu- 
tion becomes more pronounced. This apparent paradox is explained 
by the fact that when k becomes large, the higher order terms of (6) 
cannot be neglected (Table II). The resulting distributions fall to 
1/e-th their value more rapidly than the simple exponential distribu- 
tion (Figure 4), at least when z is so large that the effect of the sin- 


o--- UCb2) FOR k= 1.5 
EQ9 
— UCb2Z) FOR k= 20 


FIGURE 4 
Distribution of u(b,z) for k — 1.5 and k = 20, together with approximating 
exponentials for k = 20. : 
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gularity at z = 0 becomes unimportant. We can therefore conclude 
from this analysis that the effective longitudinal resistance of the ex- 
ternal medium is larger than the resistance calculated on the assump- 
tion that the current density is uniform. 


Stationary Distribution: Infinite External Medium 


If the kernleiter is immersed in an infinite medium, a = o, the 
calculation of the potential distribution requires the use of a certain 
generalized Fourier-Bessel integral. Weber discussed a problem of 
this sort in a later paper, (1873b), but he only treated the caseh=0. 
The generalization to the case at hand is rather immediate and in- 
volves nothing essentially new [cf. Carslaw and Jaeger, (1940) ]. 

The problem to be solved is specified as before by equations (1)- 
(5), except that now a= o. (With ®(7) chosen as an “impulse” func- 
tion, the solution will again be a Green’s function). Since the radial 
boundary is pushed out to infinity, we should expect the solution to 
involve something analogous to a Fourier integral (Watson, p. 576) ; 
in this respect the infinite problem is related to the finite problem just 
as the infinite vibrating string is related to the finite string. 

A solution of (1) satisfying the boundary conditions (2), (3), 
and at (4) ata—=ao,z—o,andr=bdis 


F(a,rje*= [hafi(a,r) +fo(a,r) Je" (27) 
where 
fila, r) = (-1)* Tolar) Yi(ab) —Ji(ab)Yo(ar)]. (28) 


From this special solution we form the general solution by super- 
position: 
w(r 2) = f B(a)F(a,r)e" da (29) 
0 
where, in order to satisfy (5), we must determine B(a) so that 


["eB@)F(a,r)da=0(0). (30) 


The expansion (29) is analogous to a Fourier integral, but is some- 
what more general. In order to determine B(a), we must find the in- 
version of (30). 

Let us multiply both sides of (30) by rF(8,7)dr and integrate 
between the limits 6b and o. Then if (r) is given by (5), we have 


[oreo Fr@,rar= at 


: zz b 


since F (6 ,b) = 2h/nb (equation 24). 


(31) 
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The left hand integral is, after interchanging the order of inte- 
gration, 


{or FG. nar ["aB(@)F(a,7)da= [°eX(a,6)B (a) da (32) 
b 0 0 


where 
X(a, B) = [or FB,1) F(a,r)dr. (33) 
b 


Since F(a, 7) satisfies the original differential equation and boundary 
conditions, we find by integrating by parts, 


X(a; fs) = 
4 
1 Pn ER a,t)a sith 
ea plimrF, ee ey aalnt ae 
If we substitute this expression into (32) we are led to consider the 
following limit: 


[Pax (a, 8) B(a) da= 
0 

(35) 
: So Bb dF (B, dF (a,7 
lim ° Se TF (4,7) a — FB 7) oe) da 
os), a BP? adr 
In this limit, the asymptotic values of F(a,7) and dF (a,7)/dr are 
required. These can be found by using the asymptotic expressions for 


the Bessel functions (Watson, p. 199) 


Fa(®) 43) cos(é —4na—i2) 
US 
4 (36) 
¥y(8) =(se) sine —bna— da) 


in equation (35). After a lengthy set of manipulations, we obtain 
from (385) and (36) 


§ bacon gh | aB(a) sin(a — B)r , 
J ao tim =| J Va p oe eet 


“a B(a) eat a B(a) sin(at Br 5 


H 
Ser caper a+ p ap “Y a+ B 


ef a B (a) K(a ) ee da 


VaBp Bs 


- (87) 
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where 
Ea) =A (CY, Ya; 65.0) iA Jae dan pane) 
Gla) =a (a Js, 8a) Aa 
Hla)=—A15 115850) ~AOG eure a 
Ki(a)= AI eo Xa, b> a) lee a 

and 


A(R, ,S8,,8,a) =Wa BR, (6b) Si(ab) — 


(39) 
h[a Ro(6 b) Si(ab) + BSo(ab) Ri (6 6)] + Ro (6 6) So(a b). 


The letters R; and S; stand for Bessel functions of either kind of or- 
der 7. 

The integrands of the last three integrals have no singularities 
on the range 0 < a < o. (The last integrand vanishes at 6 = a be- 
cause K(6) = 0.) Consequently from Du Bois-Reymond’s theorem of 
mean value (Whittaker and Watson, 1940; p. 66), the limit of each 
integral as 7 > o is zero. On the other hand, the first integral is an 
instance of Dirichlet’s integral, and so (Whittaker and Watson, p. 
176) 


f°ex(a,8) B(@)@ R= linaa| > Se 


0 ro Tt x V B oC Bb 
40 
B(S)E (8) oe 
or by the definition of F(a) [ equations (38) and (389) ], 
[Pax (a, p)B(a)da=B(6) {[h B Js(Bb) — Jo(B b)]? 
(41) 


+ Lt 8 Yas b)— Vo (8 b) 17}5 
which is the value of the integral in (32). Setting this equal to the 
expression given in (31), we have 
h-'I 


B(B) == bE(A)’ 


and if we substitute this into (29), we obtain the required Green’s 
function 


(42) 


hal 


u(r,Z) = x 
mb 
J ‘ Pages 33 
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which reduces to the function given by Weber in case h = 0. This in- 
tegral converges whenever z > 0, but is divergent for z = 0; this be- 
havior is analogous to that of the finite medium Green’s function. 
The value of u on the surface of the kernleiter sheath is found by 
putting r = 6 in equation (43). Since F(a, b) =2h/ 2b, we have 


2h2 aI 


m° b? 


u(bo,z)= x 
(44) 


oo e-a 
\j {[haJi(a b) — Jo(a b)]? + [ha Yi(a b) — Yo(a BI 


The integrals (43) and (44) cannot be evaluated exactly ; however, if 
h is very large, an asymptotic expansion of (44) can be found. Since 
(Whittaker and Watson, p. 159) 


2 0 ‘(0 "(0 

{ Be ee Os dees 5 (45) 
4 z eB 2 

the first term in the asymptotic expansion of u(b,z), in case h is 

large, is 


a 


4 
2nz’ ( 


u(b,z)~ 


while the next non-vanishing term is © 1/2° and so can be neglected. 
Unfortunately, this expansion cannot be used for small values of h, 
because (45) then becomes indeterminate. 

From (46) we may conclude that at distances sufficiently far from 
the electrodes, and for large values of h, the potential across the 
sheath falls off as 1/z. This result appears surprising in view of the 
usually accepted exponential distribution in cables. It may be argued 
that if the finite conductivity of the cable is taken into account, an 
exponential spread will result; but this is not the case, as Weber has 
shown (1873a). If the conductivity of the core is finite, the current 
spread will still be given by a sum of exponentials which reduces to a 
single exponential only when the physical and geometrical parameters 
of the system insure that the current density is uniform. 


Propagation of the Impulse 


In the present investigation we are interested primarily in two 
results: first, in the distribution of potential along the kernleiter, and 
second, in the effect of the size of the enclosing glass jacket on v , the 
velocity propagation. All of the formal conduction theories (Wein- 
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berg, 1940) start with an exponential current distribution in the deri- 
vation of the velocity formula; and while the resulting formulae dif- 
fer in their dependence on certain excitation parameters, they all de- 
pend on the length constant, 1/z , in the same way: 


Va Lee (47) 


When & is not too large, Lillie’s results seem to agree with this predic- 
tion (Figure 5) ; but if as in our case, the current distribution is only 
approximately exponential, formula (47) becomes somewhat mean- 
ingless. 

Since we are not interested in the dependence of velocity on ex- 
citation parameters, and since, in simple cases, all the theories predict 
the same dependence on external resistance, we may choose the sim- 
plest excitation equation (Blair, 1932) in calculating the velocity: 


dp/ot=u-—aop (48) 
or its steady state equivalent 
—vdp/dl=u—-op, 


where » is the excitation parameter which must exceed the threshold 
H , wis the membrane potential, ¢ = z — vt, and the voltage rheobase 
is wZ. For u we may substitute «w(b,z) from (9); if we then inte- 
grate (48) in the steady state, and set »p = H at z= 0, we find 

Patna w A; 

Oe ait Mite ot 
where the A; are defined in (7). Equation (49) is transcendental in 
v and so cannot be solved exactly, especially because of convergence 
difficulties at z= 0. If only the first term of (49) is used [cf. equa- 
tion (11)], the usual expression for v is obtained (Weinberg, 1940) 


Q) 


a (50) 


Uy — W 
Hor) 


Roughly, we would expect the velocity defined by (49) to be pro- 
portional to a “characteristic length” found by approximating the 
exact distribution (9) by an exponential over a portion of the z — 
range. Such approximating exponentials for k = 20 are included in 
Figure 4. Some ambiguity concerning the region of the u — curve to 
be approximated arises because of the singularity at z = 0; this diffi- 
culty can be partly overcome by requiring that the approximating 
curve lie above the u-curves found for k < 20 and below the u — curve 
which is the solution of the telegrapher’s equation for k = 20. 
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Experimental Implications 


An experimental study of some of the results contained in this 
paper seems to be within the range of present day techniques. For 
most of these experiments, the passive iron wire would be well suited; 
since the conductivity of its core is practically infinite compared to the 
surrounding nitric acid, our results would be more directly applicable 
to this model than to some of the collodion membrane kernleiters in 
which the core itself is an electrolyte. 


EXPERIMENTAL VELOCITY OF TRANSMISSION (LILLIE, 1925) © 
THEORETICAL [IST TERM,EQ. (Ii) | voc / — 


THEORETICAL debe DISTRIBUTION, £Q.(2/)]Vx 1. --- 
“A 


°o 2 4 6 8 10 12 14 16 16 20 22 24 


FIGURE 5 


v versus k for passive iron wire: comparison of observed velocity in passive 
iron wire with values of 1/x, and 1/x,. Note that 1/x, is too large. 


To check the general form of the potential distribution [equations 
(6) and (9)], h, the ratio of external to membrane conductivity, 
must be measured. This can be accomplished by determining, with 
probing electrodes, the potential distribution along the wire when the 
enclosing jacket just fits over the wire, and a polarizing current is 
applied through ring electrodes. Under these conditions a strictly ex- 
ponential polarization spread should be found with length constant 
given by equation (21). If the polarization is plotted logarithmically 
against distance, a straight line whose slope is the characteristic 
length will be obtained, and from this slope the value of h can be 
found, if the acid conductivity is known (cf. Bogue and Rosenberg, 
1934). The quantity h is related to the membrane structure, and so 
it should vary with the passivation time. 

For a given passivation time, h should not depend on the size of 
the enclosing jacket. Consequently the value of h found for a small 
jacket can be used to calculate several roots of the characteristic equa- 
tion (8) for large value of a; the theoretical distribution can easily 
be found from (6), and this could be checked against the experimental 
polarization. 
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A further check on the validity of the relation between velocity 
and external resistance [equation (47) ] can be obtained from the shape 
of the rising phase of the propagated action spike. It has been shown 
(Cole and Curtis, 1938; Offner, Weinberg, and Young, 1940) that the 
rising phase of the spike should be exponential with characteristic 
length, 1/x,, given by 


ty =4MCv + VV (MCv/2)? + M/s, (51) 


where C is the transverse capacitance per centimeter of the mem- 
brane. When (K) is small, (50) becomes [cf. equation 21] 


Vv — o/VM, (52) 
o being independent of jacket size. Equations (51) and (52) give 


aty/\/M = constant: (58) 


if in experiments on the variation of velocity with k , the action spikes 
are recorded, relation (53) between the moving characteristic length 
and the external resistance should be found, at least when k is small. 
When k is large, a deviation from (53) is predicted; such a deviation, 
if discovered, would tend to confirm qualitatively the theory. 

Finally, since the approximate “characteristic length” for large 
values of & is smaller than the corresponding x, (Figure 4), the ve- 
locity of propagation should fall below the curve found by Lillie (Fig- 
ure 5) when the glass jacket is large. Lillie’s results were not extend- 
ed beyond k = 5; up to this point the exact theory and the approxi- 
mate theory predict almost identical results. Data on v versus k for 
larger values of k are needed to determine the validity of the present 
considerations. 

The writer wishes to thank Professor Carl Eckart for reading 
the manuscript. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 


LITERATURE 


Bishop, G. H. 1937. “La Theorie des Circuits Locaux Permet-elle de Prévoir la 
Forme du Potentiel d’Action?” Arch. Internat. Physiol., 45, 273-297. 

Blair, H. A. 1932. “Stimulation of Nerve by Electrical Currents. I.” Jour. Gen. 
Physiol., 15, 709-782. 

Bogue, J. Y., and H. Rosenberg. 1934. “Rate of Development and Spread of 
Electrotonus.” Jour. Physiol., 82, 353-367. 

British Association Mathematical Tables, VI. 1937. Cambridge: University Press. 

Carslaw, H. S. 1906. Introduction to the Theory of Fourier’s Series and Integrals 


and the Mathematical Theory of the Conduction of Heat. New York: Mac- 
Millan Company. 


ALVIN M. WEINBERG 55 


Carslaw, H. S., and J. C. Jaeger. 1940. “Some Two-Dimensional Problems in 
Conduction of Heat with Circular Symmetry.” Proc. Lond. Math. Soc., 46, 
361-388. 

Cole, K. S. and H. J. Curtis. 1938. “Electric Impedance of Nitella During Activ- 
ity.” Jour. Gen. Physiol., 22, 37-64. 

Kalahne, A. 1907. “Ueber die Wurzeln Einiger Zylinderfunktionen und Gewis- 
ser aus Ihnen Gebildeter Gleichungen.” Zeit. f. Math. u. Phys. 54, 55-86. 
Lillie, R. S. 1925. “Factors Affecting Transmission and Recovery in the Passive 

Iron Wire Nerve Model.” Jour. Gen. Physiol., 7, 473-507. 

Lillie, R. S. 1986. “The Passive Iron Wire Model of Protoplasmic and Nervous 
Transmission and its Physiological Analogues.” Biol. Rev. 11, 181-209. 
MacMahon, J. 1894. “On the Roots of the Bessel and Certain Related Functions.” 

Ann. of Math., 9, 23-30. 

Offner, F., A. M. Weinberg and G. Young. 1940. “Nerve Conduction Theory: 
Some Mathematical Consequences of Bernstein’s Model.” Bull. Math. Bio- 
physics, 2, 89-103. 

Sasaki, S. 1914. “On the Roots of Xolit) _ eal), = 0.” Tohuko Math. 

J, (kr) J’, (ka) 
Journal, 5, 45-47. 

Smythe, W. R. 1939. Static and Dynamic Electricity. New York: McGraw-Hill. 

Watson, G. N. 1922. A Treatise on the Theory of Bessel Functions. Cambridge: 
University Press. 

Weber, H. 1873a. “Ueber die Stationaren Strémungen der Elektricitat in Cylin- 

dern.” Borchardt’s Jour. f. Math. 76, 1-20. 

Weber, H. 1878b. “Ueber eine Darstellung Willktirlicher Funktionen durch Bes- 
sel’sche Funktionen.” Math. Ann. 6, 146-161. 

Weinberg, A. M. 1940. “On the Formal Theory of Nerve Conduction.” Bull. 
Math. Biophysics, 2, 127-184. 

Whittaker, E. T., and G. N. Watson. 1940. A Course of Modern Analysis. Cam- 
bridge: University Press. 

Wilson, F. N., A. G. Macleod and P. S. Barker. 1933. The Distribution of Cur- 
rents of Action and Injury Displayed by Heart Muscle and other Excitable 
Tissues. Ann Arbor: University of Michigan Press. 


tn 


4 uy "uate 
toi a 


pes i ) AYP sf in vere 


how Le meni tied es a 


rr ye 


BULLETIN OF 
MATHEMATICAL BIOPHYSICS 
VOLUME 8, 1941 


THE DYNAMICS OF CELL CONSTRICTION DURING DIVISION 
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The approximate equation is derived for the rate of constriction 
of a dividing cell, describing the phenomenon from its early stages. 
The equation previously derived by G. Young for the case when the 
constriction has already considerably progressed is obtained as a lim- 
iting case. a ives 


The mechanism of cell constriction during division has been 
studied mathematically by G. Young (1939; c.f. also Rashevsky, 
1940). An equation for the rate of constriction has been derived for 
the case when the constriction has already proceeded rather far. This 
equation has proved to agree rather well with available experimental 
data (Rashevsky, 1940, chap. iii). For the early stages of constric- 
tion no mathematical treatment has yet been available. By using the 
exact solution of the diffusion problem for an ellipsoid of revolution. 
and applying Betti’s equation to different regions of the ellipsoid, G. 
Young (1939) was able to show that a constriction will set in, but 
no expressions for its variation with respect to time was derived for 
those early stages. 

It is the purpose of the present paper to derive a crudely approxi- 
mate, but general expression, which would describe the process of 
constriction in all its entirety. We shall follow the method used by 
G. Young for the case of already advanced constriction. One addi- 
tional factor will be introduced, by way of correction to Young’s ar- 
gument. 


be 


FIGURE 1 


Consider the situation represented in Figure 1, in which the no- 
tations are obvious. This corresponds to G. Young’s case (Young, 
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1939; Rashevsky, 1940, chap. iii). We shall now drop the restriction, 
used previously, namely that the radius 7 of the neck is small, com- 
pared to 7, and r,. Before the constriction sets in, r may be as large 
as the smaller radius 7, of the cell. 

We have as before the force 


(1) 


which pulls on the neck. Even when the neck is very narrow, equation 
(1) is only a very rough approximation, obtained by neglecting the 
disturbance of the diffusion field of one sphere by the other. The same 
degree of approximation equation (1) may be used for the case, in 
which constriction has not yet set in. We may consider that the re- 
pelling force of each half of the cell is roughly the same as that of 
two closely adjacent spheres, with radii 7, , where 27, is the average 
over all width of the cell. 
To the force F is to be added the force 


HART vy (2) 


due to the lateral surface tension of the neck. In addition however a 
force is acting on the “ends” of the neck, namely the force 


2 
P=-2" xar (3) 
Te 


due to the capillary pressure in the “spherical” parts. This term has 
not been taken into account by G. Young. 
Furthermore there is a pull upon the “ends” of the neck equal to 


2nrycos@=2aryX— 7? (4) 


due to the tangential component of the surface tension of the “spheres” 
at the edge of the neck, where the latter is connected to the “spheres”. 
This force may be best combined with the force (2), giving as the 
total force due to tangential surface tension. and directed along the 


axis of the neck \ 


~2ary(l-—). (5) 


When the neck is narrow, r < 7., we have only the term used by G. 
Young. When on the contrary, the constriction has not yet begun, 
this term is zero. 


Altogether the forces (1), (3) and (5) are equivalent to an out- 
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ward directed pressure p,, applied to the “ends” of the neck, and 
equal to 


Fr." — 2nyr? — 2ayr(r2" —7r)  F — Quyr 


oS i, eT od 
At the “sides” of the neck we have the capillary pressure 
hat? 
Peo r ’ (7) 


as before. Applying now Betti’s formula to the neck, and neglecting 
again the volume forces, which will be small due to the small length / 
of the neck, we find, just as before (G. Young, 1939; Rashevsky, 
1940) 


ar 


y deoel 
ZieaeG a a 6 rs ey 
This is identical with Young’s equation. However, by the nature of 
its derivation, it should hold even for the early stages of constriction. 
In fact, the comparison with experimental data (Rashevsky, 1940, 
chap. iii), which gave a rather good agreement, was carried out over 
a much wider range of values of 7, than would have been permitted 
by the restriction to very advanced constriction. 

The essential difference between our derivation and Young’s is 
that we do first correct Young’s argument for the capillary pressure 
term (3), which correction should hold for any degree of constriction ; 
and second, introduce the term (4), which plays a role only when the 
constriction is small. The two terms compensate each other. If we do 
not introduce correction (3), but only correction (4), we would find 

Bee (nae AM Rite ty, (9) 
dt. .6y 377° 6Gayr 
which would reduce to Young’s equation only for very small values 
of r. 

One feels somewhat uncomfortable however in introducing a 
sharp angle between the “neck” and the “spheres”. And if we con- 
sider that actually such an angle does not exist, but that the meridion- 
al line of the cell is a smooth curve, the whole argument which leads 
to expression (5) becomes rather doubtful, and even the original ex- 
pression (2) used by Young loses its meaning. We shall therefore 
treat the problem in a different way, remembering that a positive 
pressure at the “ends” of the neck is equivalent to a corresponding 
negative pressure at the “sides”. 
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We shall consider, as is actually the case, that the neck is nega- 
tively curved in the meridional plane, with an average radius of cur- 
vature?”. Then we have for the pressure p, at the “ends” of the neck: 


eect 2 


a ee 10 
ue gee oe 
and at the “‘sides”’: 
Y v 
ee 11 
De al x (11) 
Since, from Betti’s formula 
ldr 1 
fs ee eRe — pz), 12 
ae 67 (91 — Pz) (12) 
therefore 
Taek F 2 T; 
ast ar ental fo 


The radius 7’ may be approximately calculated in the following way: 
The quantities 7,” and r.” are approximately connected with the 
overall length of 27, and the overall width 27, of the cell by 


Sh fo a (14) 
while the half length / of the neck is given by 
‘= ae > a Vo =T7; ia 27> ° (15) 


Expression (15) holds only for 7, > 27,. A more general but crude 
approximation may be made in the way indicated in Figure 3. 
When the neck is very wide, the notion of length of the neck be- 
comes somewhat indefinite. We shall use equation (15) as defining 
the length of 2/ of the neck for any value of r. 
For r = r, the radius 7” is infinite. This at first seems to contra- 
dict the assumption made in the theory of cell elongation, according 


FIGURE 2 
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to which the mean curvature at the sides of the cell is taken to be 
1/r. + 1/71, which means rv’ = — r,. We must remember however, 
that in the latter case we refer to the average mean curvature of the 
“sides”. At present however we consider only the curvature in the 
meridional plane near the equator. 

If during the process of constriction 7, and 7, remain approxi- 
mately constant, then 7’ decreases from © to a value approximately 
equal to 1 (Figure 2). This value is roughly reached when 7, — 7 
becomes equal tol. Between r. — r = 0 and 7, — r=1, the value of 
r’ may be calculated from Figure 8, using the relation 7’/y = y/2z, 
and is found to be equal to 


, 


1 (rm. —r)? +P 
(= SS 


% Ti iF (16) 


From r = 7, to r= 72 — I the value (16) should be used in equation 
(183). For 7 < r. — 1, the value 7’ = 1 is to be used. For very small 
values of 7 the first term of the right hand side of (15) prevails, and 
equation (15) reduces to that of Young. The condition for the con- 


FIGURE 3 


striction to occur at all is obtained by requiring dr/dt < 0 for r= 712. 
This gives 


Rey ats, (17) 
or, because of (1) and (14) 
RT uwaqr? 
6M jay (Tr rm 12)? 


a relation that can be satisfied with plausible values of the constants. 
If relation (18) is satisfied, then the constriction once started 
for r=7., will proceed until complete separation. 


Sle (18) 
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The author is indebted to Mr. H. D. Landahl for checking the 
calculations. 
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A THEORY OF STEADY-STATE ACTIVITY IN 
NERVE-FIBER NETWORKS: 
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As an essay towards the determination of the effect of structural 
relations among nerve fibers upon the character of their activity, 
preliminary consideration is given to the steady-state activity of some 
simple neural structures. It is assumed as a first approximation that 
while acted upon by a constant stimulus, each fiber reaches a steady- 
state activity whose intensity is a linear function of the applied 
stimulus. It is shown by way of example that for a simple two-fiber 
circuit of inhibitory neurons knowledge of the stimuli applied to the 
separate fibers does not necessarily suffice to determine uniquely the 
activity that will result. On the other hand, there are deduced cer- 
tain restrictions on the possible types of activity that may be consis- 
tent with a given pattern of applied stimulation. 


Studies in the mathematical biophysics of the nervous system by 
Rashevsky (1938, 1940; also various papers in Psychometrika and the 
Bulletin of Mathematical Biophysics), Landahl, and Householder (in 
the same periodicals) have presupposed certain equations as describ- 
ing the dynamics of individual nerve fibers, and have deduced from 
these presuppositions the equations which would describe the behav- 
ior of certain nerve-fiber complexes of specified structure. These equa- 
tions express the time course of the response of the efferent fibers of 
the complex in terms of the intensity of the stimulus applied to the 
afferent fibers, so that if the response of the organism depends in a 
known fashion upon the response of these efferent fibers, and if the 
intensity of the stimulus which acts upon the afferent fibers is given 
in terms of the external stimulus, then the response of the organism 
can be said to have been derived as a function of the external stimulus. 

Hitherto all derivations have been made on the basis of special 
structural assumptions with no attempt at greater generality. Never- 
theless, if we neglect conduction time (along the fiber), which is cer- 
tainly legitimate to at least a first approximation, the behavior of any 
complex of nerve fibers must depend alone upon the dynamic proper- 
ties of the individual fibers and the structural relations among these 
fibers. (The possibly variable physiological state of the organism as a 
whole, which may well effect the activity of some or all of the constit- 
uent fibers, is here assumed to constitute a part of the acting stimu- 
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lus.) Hence there must necessarily exist certain general propositions 
which express the activity of any nerve-fiber complex in terms of the 
dynamics of the individual fibers and of the structural constants. The 
present paper serves as an initial step towards the search for such 
propositions. 

Naturally steady-state activity under constant stimulation is the 
simplest to consider; with this we therefore begin and to this we con- 
fine the present paper. That is, if a constant stimulus is applied for 
a period of time to each of certain fibers of the complex, whatever the 
fibers may do in the first few milliseconds, or even seconds, we sup- 
pose they settle down shortly to a steady response which persists un- 
changed as long as the stimuli are left unchanged. It is this steady- 
state behavior that we investigate here. 

We assume, as a first approximation, that the steady-state activ- 
ity of each fiber is a linear function of the intensity of the stimulus 
applied to it. Since this necessarily vanishes when the stimulus falls 
below a certain intensity which we here call the threshold (and which 
might conceivably exceed the true threshold), it follows that the re- 
sponse is to be taken as proportional to the amount by which the 
stimulus exceeds the threshold. 

If a nerve fiber is stimulated at several points, and if the steady- 
state response is additive, as we assume, it is no restriction to think 
of the single fiber as replaced by several, one passing from each den- 
dritic to each axonic synapse of the fiber. We therefore suppose that 
each fiber of the complex has but a single initial or afferent or den- 
dritic synapse, and but a single terminal or efferent or axonic synapse. 
We further suppose that given any two synapses, say s, and s., not 
more than one fiber has s, as its initial and also s. as its terminal 
synapse, and moreover that no fiber initiates and terminates in the 
same synapse. This is no real restriction, for any single fiber can be 
replaced by two or more in a chain, and the assumption is made for 
notational convenience only. 

For if there are n synapses, designated by si, 82, +++, S,, and if 
there is a fiber leading from s; to s;, we may denote it by [77], since 
there can be no more than one. Then at any given synapse, s,, 
suppose there is applied a stimulus of constant intensity, S; , positive, 
negative (inhibitory), or zero, from a source external to the complex 
itself and therefore independent of its activity. If constant stimuli 
are being applied at the same time to other synapses of the complex, 
and if there are fibers [17] terminating in s, , then some of these may 
have been set into activity. If so the effective stimulus acting upon 
any fiber [7 k] initiating at s; will be, not S;, but the sum of S ; and 
the amounts of « — 7 produced by these fibers. Let n; represent the 
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value of this sum after the system reaches its steady state. In the 
final steady state for the system each fiber [7 7] will be acted upon by 
an effective stimulus of amount 7; and produces e — j in an amount 
equal to ai;(7i — hi;). Here h;; is the threshold of the fiber [77], and 
therefore essentially positive, while 


0 when 9: < hi;, 
O55 = (1) 
Qi; when Hi > hi; - 


The parameter a;; is a characteristic of the fiber which may be called 
its activity parameter. It is assumed to be different from zero, but it 
may be either positive (when the fiber is excitatory in character) or 
negative (when the fiber is inhibitory in character). Thus, since at 
each synapse the value of is equal to a sum of the S applied from 
without the system, and the amounts of « — j produced by the active 
afferent fibers, it follows that 


03 = S; + = ais (ni — his). (2) 


The summation is extended over all values of 7 for which there exists 
a fiber [77], and the coefficients a;; are defined by (1). Since the a’s 
are themselves functions of the 7’s, these equations do not constitute 
a simple linear system, but rather present, for any given set of values 
of the S;, a set of linear systems, one for each choice of a set of a’s 
to be taken as non-null, and among this set at least one system will 
presumably be such that the solutions 7 will be consistent with the 
conditions (1). But there is no assurance in advance that only one 
such choice of the non-null a’s will yield 7’s as defined by (2) which 
are consistent with (1). In fact, a simple example will be exhibited 
below in which there are certainly more than one—at least for suit- 
able S — values. In other words, the S — values may not define 
uniquely the steady state activity of a given nerve-fiber network. In- 
stead, it may be necessary to know in addition the order in which the 
stimuli S; are applied, or the rapidity with which the various fibers 
of the complex are capable of reaching their own steady states. 

For a given system of nerve fibers let the values S; of the stimuli 
externally applied be called the stimulus pattern SP. Let the values 
of the 7; be called the excitation pattern HP. And let the values of 
the a;; constitute the activity pattern AP. Hence the general problem 
is to determine, for any given SP, every possible AP (or EP) that is 
consistent with (1). We remark in passing that when the structure 
of the system and the values of the parameters aj; and hj; are known, 
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then for any EP, the AP can be determined by (1) and the necessary 
SP by (2), both uniquely. 

We note now that for a simple chain of fibers [12], [23], ---, 
[n—1,n] the AP and the EP are uniquely and easily determined by 
the SP. For the activity of [12] and hence the value of 7. depends 
upon S, alone; the activity of [23] depends upon 72 alone; and so on 
to y7,. We therefore no longer need to consider a simple chain. If the 
system contains a simple afferent chain leading, say, to s;, we may 
disregard it and consider S; as including the contribution of this 
chain. If the system contains a simple efferent chain leading, say, 
from s;, we may disregard it since the activity or non-activity of the 
fibers of this chain can have no effect on the rest of the system. Hence, 
we are led to consider only such systems containing synapses 
8;(7 =1,---, ”) such that every synapse joins at least one afferent 
fiber [77] belonging to the system to at least one efferent fiber [7 k] 
belonging to the system. We suppose, moreover, that given any two 
synapses s; and s; there is at least one chain [791] , [71 Jo] , -:: [7 k] 
leading from s; to s, and also one leading from s, to s;. Hence, the 
system does not consist of two separated subsystems, nor does it con- 
sist of two subsystems, say A and B , of such a nature that while there 
are fibers or chains leading from synapses in A to synapses in B, 
there are none leading from synapses in B to synapses in A. Sucha 
system, whose every synapse joins an afferent to an efferent fiber, 
and whose every pair of synapses is joined by chains going both ways, 
we shall call a fiber network, and these alone we consider hereafter. 
The discussion of any other system can be reduced to the discussion 
of one or a chain of such networks. 

The simplest possible network is a simple circuit of fibers [12], 
[23], --- , [n,1]. Discussions have already been given of the non- 
steady-state activity of a single self-exciting fiber (Landahl and 
Householder, 1939), and of a pair of mutually exciting fibers (House- 
holder, 1938), both excitatory. It is easy to see that 

Tf all fibers of a simple. circuit are excitatory, then not more than 
one AP (EP) can be consistent with any given SP . 

For a simple circuit we may simplify the notation somewhat, 
dropping the second subscript on the symbols a, a and h » SO that 


Qj = Dijit » Qi = Aiin , h,= Nijivs : (3) 


(All subscripts may be reduced modulo n.) Then define the quanti- 
ties 


Ga =S; vee 9 


o™ =e + A 0"), (4) 
7 % 4-1 c 
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Evidently 
of? oe ate Qji-1 °+ + Di-a oy aie (5) 


Also let 
Y¥i=n—h. (6) 
Then for the simple circuit the equations (2) reduce to 
Yi oi + in Yin. (7) 


If the fibers are all excitatory, the activity of ‘any one cannot de- 
crease the activity of the succeeding one. Hence if every o; < 0, no 
fiber will be active. If every a; > 0 then every fiber must be active, 
every a; — a; , and (7) can be solved by ordinary methods. In the re- 
maining alternative there will be some positive and some negative o’s. 
Suppose o; < 0. If [7,7 +1] becomes active it can be only after and 
because [7 — 1 , 7] has become active. Now starting with each oj 
which follows a non-positive o; , form the sequence 


(1) (2) 


Oj = ofr» Tj12™ 2+, Tjxv™ y (8) 


terminating with the first non-positive term (which may be o;,; it- 
self), or else with o;™. If, among these sequences, there is one which 
contains all n terms, then either the last term, o;, is non-positive, in 
which case the solutions y of (7) are given by the members of this 
sequence, 

Yjra = oj10™, (9) 
or else the last term a; is positive and all fibers are active, the HP 
being given by (7) with a;—a;. But if no sequence has nv terms, then 
eliminate every sequence whose initial member o;,, has a subscript 
7 + 1 which is the same as the subscript of some member of another 
sequence. There will remain exactly n of the o’s, 


oa, 6) 4 oo (%) pitty on™ 


whose superscripts have the property 


Qin — 0; +1 when oi >0, 


10 
Cag — 1 when o;‘%) <0. eM, 


Then the only AP consistent with the given SP will be given by 
Yi oi'*, (11) 


Return now to the case when all fibers are excited, which occurs 
whenever a sequence (8) exists having no negative members. The 
determinant of equations (7) in the case when all a’s are non-null is 
equal to 
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A= 1 — AyAy +--+ On , (12) 
as is easily verified, and the solutions in this case are 
Yi=oai™/A. (13) 


But every y; > 0 and at least ojin > 0. Hence 4 > 0, which means 
that the parameters a; cannot be too large or the fibers too active. 
This apparent paradox simply means that if the activity parameters 
a; are too large and the fibers too active, then if all fibers were set 
into activity at once the excitation would build up beyond the point 
at which the linear approximation can be retained. Hence 

For a circuit of only excitatory fibers the linear approximation 
does not give a steady state for the case when all the fibers are active 
if the product of the activity parameters a; exceeds unity. 

A further lemma is almost immediate and will prove useful in 
later investigations. 

If, for a given SP applied to the synapse of a simple circuit, there 
are two distinct AP’s both consistent with this same SP, then every 
a, which is zero in the one must of necessity be non-zero in the other. 
That is to say, if a given set of constant stimuli S; applied to the 
synapses of the circuit is capable of leading to either of two distinct 
AP’s according to the order or manner of application of the S;, then 
every fiber which is inactive in the one case is necessarily active in 
the other. 

This is vacuously true when the circuit contains no inhibitory 
fibers. But in any case if there is an inactive fiber in a given AP, then 
its removal from the circuit would have no effect upon the AP or the 
EP of the remaining system, which is now a simple chain. But for a 
chain the SP determines uniquely the AP and the EP, so that the 
statement of the lemma follows. 

We conclude these preliminaries with a brief consideration of the 
circuit of two fibers. We have to consider the quantities 


(1) 


oe Shee 
O71 Nit g o,) =o, + A, oo, 


oe") —= 02, a2") a OD a Oy 01. 


If o, and o, are both non-positive no activity is possible, but if at least 
one is positive then at least one fiber must be active. Hence on account 
of the lemma just proved, if for a given SP there are two distinct AP’s 
both consistent with this SP, they can only be described by the values 
(a,,0) and (0, a) respectively of (a, , a2), for it can be shown that 
the AP(a,, a2) is not consistent with either of these. The first AP 
is associated with the EP y, = 0, > 0, y2 = o2 + a; «1, < 0; the second 


with Y2= 02 > 0, Y1 = 0; + G op < 0. We inquire, therefore, whether 
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for any choice of the S; , or, what comes to the same thing, of the o;, 
the four inequalities 


q- 0; re i ats ee Oe 
a0, Ga 8s oy = 0, 


may hold simultaneously. If a, <0, a. <0, anda, a > 1, then in the 
(o;, o2) — plane it is readily verified that the ray o. + a, o, = 0 in 
the first quadrant lies above the ray o, + ad o, = 0 in the same quad- 
rant, and that any choice of the o’s in the angular region between 
these rays will satisfy the above inequalities. For such a circuit so 
stimulated either of two different activity patterns (a,,0) and (0, a) 
may occur. But if either of the three inequalities on the a’s fails then 
the SP determines uniquely the AP and the EP. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of the University of Chicago. 
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A special case of a problem discussed in a previous paper is 
treated in greater detail. An equation in the three variables, errors, 
trials and number of possible choices, is developed and compared with 
the results of an experiment performed under conditions closely ap- 
proximating those required for the development of the equation. The 
agreement is excellent. 


The development of a mechanism for discrimination and condi- 
tioning has been given elsewhere, and relations have been derived 
for special cases (Landahl, 1941, hereafter referred to as A). Here, 
we wish to consider a special case, involving only the variables, trials, 
errors, and the number of possible choices, N. The variation with N 
was discussed very briefly in A. We continue to use the same notation. 

Instead of a choice between two stimuli, let there be a choice from 
among WN stimuli, --- S’;---. The stimuli are assumed to have the same 
initial stimulus values, so that experimentally all bias is averaged out. 
The general case can be treated but with considerably more difficulty. 
Let ¢ represent the initial excitation factor value of any one of the 
stimuli acting on the mechanism. Let one of these stimuli S’, be posi- 
tively associated with a stimulus S,, outside the primed group of 
stimuli, so that if c is the number of correct responses, then ¢,, the 
stimulus value for S’., is given by [A, equation (11) ] 


&=& + be, (1) 


where b = B f(t) in equation (11) of A. 

Let us next assume that the stimuli of the primed group, among 
which a choice is made, other than S’., are negatively conditioned to 
S.. If #’ is a positive constant analogous to b , and if w; is the number 
of responses to some item S’; produced by S., then the stimulus for 
S’; is given by 

Ew, — €o — B Wie S f (2) 
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We wish, however, to consider only the correct and wrong re- 
sponses in general and not wrong responses to particular stimuli. 
Hence, we obtain the average of the stimulus value for the wrong 
stimuli, and designate this value by «,. Summing equation (2) over a 
for the N—1 stimuli, we have for the average value e , 


1 B 
== a a= _- ie 3 
Ew Noi: Eo woi2” ( ) 


As, © w; = w, the total number of wrong responses, equation (3) be- 
comes 
Ey = & — Bw, (4) 
where $8 = f’/(N—1). Here f is not analogous to b of equation (1). 
When we neglect failures to respond, we have 
N= Ow’, (5) 


where n is the number of trials. 

If P., is the probability corrected for chance that the choice is 
correct, we have from previous development [Landahl, 1938, equa- 
tions (10) and (15) ] 


log (1 — Pa) =—K(ée — vw) =—k(bet+ pw). (6) 


The relation between P., and P,., the probability for a correct re- 
sponse (not corrected for chance) is given by [Landahl, 1939, equa- 
tion (29)] 


we WN i 
Pa= ay (Pe w: (7) 
But we also can write instead of P.., 
dc 
Cm adn 7 (8) 


since we identify P,. with the instantaneous rate of making correct 
responses. Thus, introducing (8) into (7) and the result into (6) 
and using (5), we obtain 

dw 1 dc N-1 


ae 7 = W e-k(Bw-bw+bn) | (9) 


Integrating equation (9) and setting the initial conditions as 
¢ = w=n= 0, and solving for w, we obtain 


1 Nb 


Oe B) = WAT O= Hers N= 1ETE (19) 
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We shall consider two special cases of equation (10). First let 
us neglect the effect of the negative conditioning so that b=, =0. 
Then 


N 


1 
SA Og et eee 8 = 
mee iiica(Nesleru oo 


(11) 

Under certain conditions one might plausibly expect that f=bdb, 
that is, the positive and negative learning are equally effective. In 
this case 6’ = b = (N—1)8. Then 


peggy Los N es 
Ne 2b (N — 2)e¥n +92 Sea eas, 


The equations which have been used hold well enough for prob- 
abilities not much greater than chance. For large probabilities of a 
correct response, error is introduced by assuming that we may use 
&c — €w. This value is too large as only those values of ¢,, near e, have 
much influence. Thus P. is overestimated. On the other hand, for 
large P, , the approximation of a Gaussian distribution by an absolute 
value exponential curve [Landahl, 1938, equations (6) and (7)] un- 
derestimates the value of P.. As these two errors, which occur at 
large values of P., tend to cancel out, we may hope for a rather good 
Tt. 

We shall now discuss the relation of the above equations to ex- 
periment, but first mention reasons for using the particular proce- 
dure which is described below. One is that a solution was found for 
this case, and the experiment was thus adapted to fit as closely as pos- 
sible the theoretical assumptions. But also, as the development in 
A is made with particular reference to comparative data, and as we 
are here treating a special case of the more general development we 
have recognition rather than recall learning. The procedure could 
perhaps have been adapted to obtain both results at once, as also the 
theory; but this was not done. 

One hundred and twenty pairs of four letter words were pre- 
pared in groups of four, eight and twelve—eight groups of four, five 
groups of eight, and four groups of twelve. The words were paired 
with some care, and there were no duplicates in the two hundred and 
forty words. Of each pair of words, we may refer to one as the stimu- 
lus word (S) and the other as the response word (S’). For each of 
the groups of four, eight and twelve response words two cards were 
prepared, each containing the same (four, eight or twelve) words, 
arranged in a different order. Alternate cards were used so that the 
effect of position on the response card would be largely eliminated 


(12) 
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and also to prevent the subject from using some system in learning 
based on position. The stimulus words in each group were written 
on individual cards. 

One subject was given six sittings at each of which twenty pairs 
of stimuli were presented sixteen times. Only one sitting was given 
in one day. Since the rate of presentation of items was approximately 
constant, the average time between an item and its repetition was the 
same for each of the groups (four, eight and twelve), and the same 
in each of the sixteen trials (about seven minutes). 

For a particular sitting the subject was given a group of eight 
and a group of twelve pairs of words. For any trial the eight or 
twelve group was chosen at random to be first. If the eight group 
was chosen first for a particular trial, one of the two response cards 
(identical except for arrangement) was presented to the subject. A 
card containing a stimulus word frem the corresponding group was 
chosen at random from among the eight stimulus words and present- 
ed. The subject had been-told that for each stimulus word there was 
one word on the response card which was to be associated with it, 
and that he was to determine the correct response word in each case 
by trial and error. He was given but one chance to respond to a 
stimulus word in a particular trial, each of which consisted of twenty 
responses. He was to remember to give correct responses and to 
avoid wrong responses on successive trials. The response was ver- 
bal. The experimenter said “right” if the response word was paired 
with the stimulus word on the answer sheet for that group, if not, 
he said “wrong”. After about ten seconds the next stimulus was pre- 
sented. Immediately after the eighth item had been presented, the 
twelve group was similarly presented. One trial was thus completed 
in about five minutes. As the time the subject took to respond varied 
considerably with trials, the period between trials was varied, largely 
compensating for the former. This period was about two minutes, 
and was spent in reading or other activity of interest to the subject. 
The procedure was carried out to sixteen sets of trials of which the 
last three trials were errorless trials, except in two cases. The two 
exceptions were given extra trials. 

The constants b and f’ are functions of the time between repeti- 
tions as well as the average number of intervening presentations (in 
this case twenty). But this was constant throughout the three groups. 
However, these constants may also depend on the number in the 
group—four, eight or twelve. This would be possible, if within each 
group there was some sort of mutual inhibitory effect. The effect 
might also be present as an influence of one group upon another in 
the same sitting but the latter would be averaged out due to the vari- 
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ous combinations of three groups and to the random order of presen- 
tation. But the effect of stimuli within the group upon one another 
can show up if present. Hence b = b(N — 1) in general. Now, if we 
assume some mutual inhibitory effect (cf. Rashevsky, 1938, Figure 
54), 6 will decrease with N. But b must be positive or zero. Let us 
therefore somewhat arbitrarily assume the relation 


wy 
b= 7° tN, (13) 
It must be kept in mind that the particular form of equation (13) is 
not as yet derived from the mechanism, but its general shape can be 
deduced. We shall find that we could almost as well use an equation 
of the form b = 7/(1 + CN)? or b = 7” — 6" N. Equation (13) is 
somewhat more convenient and fits the data slightly better. However, 
the accuracy of the latter, and the presence of only three values of N 
do not warrant any choice as to the form of the b(N) relation as long 
as it is a decreasing function. With more accurate data containing 
more values of N , this relation could be more accurately obtained. 
As the subject was given an equal amount of time to remember 
a correct or wrong response, one might expect that b = f’, so that 
equation (12) should fit the data. However, the subject’s report that 
he did not attribute much value to remembering the wrong responses, 
although instructed to do so, would justify the use of equation (11). 
It was found that equation (12) did not fit any one of the curves very 
well, so that equation (11) was used. In order to fit the three curves 
simultaneously, we must also use equation (13). Introducing this 
equation into equation (11), we have a two parametric surface in the 
variables w,n, and N: 
W = ee log ye, = 9 (14) 
7 (N=-A)je™m +1 


the parameters being 7 and C. 

In Figure 1 are graphed the curves w(n) as given by equation 
(14) for N=4,N=8, and N = 12, using the values y = 1.15 and 
¢t =0.098. The crosses, open circles and solid circles are the data for 
the respective groups, four, eight, and twelve. Each cross is a value 
based on thirty-two responses; each open circle is based’on forty re- 
sponses; and each solid circle is based on forty-eight responses. The 
ordinate represents the cumulative wrong responses for one stimulus ; 
the abscissa represents the number of the trial; and N is the num- 
ber of possible responses in the group. The agreement is better than 
would be expected. Even much poorer agreement would be meaning- 
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ful, since, we are treating a special case of a more general situation. 
A test, using different material, under different conditions, especially 
comparative data, would be desirable. Only then, if sufficiently close 
agreement is found, could the development be considered as satis- 
factory. 

By taking a fixed value of N , we may find the relation between 
bk and such variables as time between presentations and effect of re- 
ward; and by using the more general equation (10) the ratio of b to 
6’ can be obtained for varying degrees of reward and punishment. 

It may be noted that it is the value b& which enters as a para- 
meter—the product of b, the coefficient of conditioning, and k, the 
parameter of the distribution curve, k being approximately the re- 
ciprocal of the standard deviation. Thus, although b is independent 
of discrimination and k is independent of conditioning, we find that 
the final rate of learning is determined by the product of these two 
parameters, and b itself is a product of two parameters. Such a situa- 
tion occurs frequently as in the various cases discussed in A. 

If such data are available from a number of experimental situa- 
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tions, it would be necessary to analyze the values involving products 
to obtain the separate parameters. If the data are available for a 
group of subjects, one can express the separate parameters as well 
as the products as deviates from the means. Expanding the result 
and retaining only the terms of first degree in the deviates, the ex- 
pression for an experimental value is given as the sum of terms each 
of which is a coefficient times the deviate value of a basic parameter. 
This may be summed over all the parameters as the deviates are zero 
for those which are absent. Generally one would expect a particular 
experimentally derived value to involve only a limited number of the 
basic parameters, as is the case above. Thus we have the basic equa- 
tion of factor analysis as well as a requirement of simple structure 
(Thurstone, 1935). The factor method might well be useful if such a 
stage is reached and if a simpler method is not possible. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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SCOPE OF THE BULLETIN 

1. The Bulletin is devoted to publications of research in Mathe- 
matical Biophysics, as contributing to the physicomathematical foun- 
dations of biology in their most general scope. 

2. Papers published in the Bulletin cover physicomathematical 
theories as well as any other mathematical treatments of biological 
phenomena, with the exception of purely statistical studies. 

3. Mathematical studies in physics or in borderline fields in 
which a direct connection with biological problems is pointed out are 
also accepted. 

4. Emphasis is put upon the mathematical developments, but a 
description and discussion of experimental work falls also within the 
scope of the Bulletin provided that description or discussion is made 
in close connection with mathematical developments contained in the 
same paper. 

5. Outside of the scope of the journal are papers of purely sta- 
tistical nature or papers concerned only with empirical equations. 


PREPARATION OF MANUSCRIPTS 


All manuscripts should be typewritten double space. Equations 
should be numbered on the right consecutively. Do not use such nota- 
tions as “equation 2a” or, “equation 5’” etc. References should be all 
given at the end of the paper, arranged in alphabetic order by authors 
and for each author chronologically, following exactly the style used in 
the Bulletin. In text, reference should be made by giving in paren- 
theses the name of the author followed by the year of publication. In 
case of several publications by the same author in the same year, use 
notations “1940a”, “1940b”, etc. 

In writing equations, slanted lines should be used wherever pos- 
sible. 

Every paper is to be preceded by a short abstract. 

Drawings should be prepared in a professional manner on white 
paper or tracing cloth following as closely as possible the style used 
in the Bulletin. They should be drawn approximately twice the scale 
of the finished reproduction. Lettering should be made by means of a 
lettering guide. If desired by the authors, drawings may be made at 
the Editorial Office according to author’s instructions, at cost. 
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